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Abstract 

Discrete Fourier transform (DFT) is the base of modern signal or information processing. 1- 
Dimensional fast Fourier transform (ID FFT) and 2D FFT have time complexity O(NlogN) and 
0(N 2 logN) respectively. Quantum ID and 2D DFT algorithms with classical output (ID QDFT 
and 2D QDFT) are presented in this paper. And quantum algorithm for convolution estimation 
is also presented in this paper. Compared with FFT, QDFT has two advantages at least. One of 
advantages is that ID and 2D QDFT has time complexity 0(y/~N) and 0(N) respectively. The 
other advantage is that QDFT can process very long signal sequence at a time. QDFT and quantum 
convolution demonstrate that quantum signal processing with classical output is possible. 
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I. INTRODUCTION 



A. Introduction of Discrete Fourier transform (DFT) 

DFT is the base of modern signal and information processing. No DFT, no modern signal 
and information processing 
Let 
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Wn defined in EqJT]is called Fourier transform matrix 
1-Dimensional DFT (ID DFT) [1] is defined as 



' W »x ^ 



Wn-i is the line vector of the matrix. The matrix 



Wi • x 



W N _i • x 



(2) 



, where vector x = (xo, x±, x^ y • • • ,xn-i) t , c = (co, ci, C2, • • • , ctv_i) t (T denotes the 
transpose of vector), and "•" denotes the inner product between two vectors. Each compo- 
nent of vector c is called Fourier coefficient. 



For any input matrix F = [f. 
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2D DFT 



in 



is defined as 



C = W N FW. 



N 



(3) 



, where C denotes the matrix of DFT coefficients. 

DFT has two important properties, that will be applied to design quantum algorithm in 
this paper. One property is that DFT is energy conservation transform (i.e., | \x\ | 2 = 1 1 c 1 1 2 , 
or ^2(xi) 2 = |cj| 2 ). The other property is that, typical data sequence, such as digital 
image, has high redundance, many of its DFT coefficients have values close to zero, and 
these coefficients can be discarded without seriously affecting the estimated value of ||af|| 2 



(e.g., without seriously affecting the quality of the restored image). Therefore, inverse DFT 
acting on the fe' pproximatively. 
FigfTJ shows the ] 




FIG. 1: The Illustration of the Two Properties of DFT that Applied to Design Quantum 
Algorithm in This Paper: fig. A denotes the original data of the top-left image block with size 
4x4 of image Lena fl^ |. Fig.B denotes the 2D DFT coefficients of the image block. DFT has two 
properties. One property is that DFT is energy conservation transform. The other property is that, 
only few DFT coefficients is not close to zero, and the inverse DFT acting on these coefficients can 
restore the original data approximatively. The two properties are applied to design quantum DFT 
algorithm in this paper. 

Fast Fourier transforms (FFT) with time complexity O(NlogN) and 0(N 2 logN) for ID 
and 2D DFT respectively were presented at 1965 [1]. And the desire for designing more fast 
DFT algorithm is still very strong. However, there is no more fast algorithm is presented 
up till now. Classical computation can not improve the efficiency of DFT any more maybe, 
and the principle of classical computation becomes the bottleneck of designing fast DFT 
algorithm maybe. Can quantum computation be applied to design more fast DFT algorithm 
with classical output to break the bottleneck? 



B. Introduction of Quantum Computation 

Introduction of Quantum Fourier Transform (QFT): QFT J2I, 3] on an orthonormal 

N-l 

basis |0), |1), • • ■ , \N - 1) is defined to be a linear operator QFT\j) = -L e 2lTijk/N \k). 
QFT is the key of the famous Shor's order- finding and factoring algorithm [2] . QFT has time 
complexity 0(log 2 N). However, QFT is not suitable to signal and information processing 
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t of DFT denned in Eqj2] (or Eqj3]) can not be generated and measured out 



because the resu 
by this QFT t 

Introduction of Grover's Algorithm: Grover's algorithm [6| solves the problem of 
searching for an element with a unique index i in a list of N unsorted elements, similar 
to searching a database like a telephone directory when we know the number but not the 
person's name (jj. Grover's algorithm has time complexity O(VN) 3|. Long proposes a 
modified Grover's algorithm, that has the probability of success 100% even for the case 



that the number of elements is very small 



7|. Boyer, Brassard, H0yer, and Tap present 



the modified Grover's algorithm named BBHT algorithm in this paper for the case that the 
number of solutions is unknown [8|. BBHT algorithm is a very smart algorithm because it 
saves many quantum circuits. 

Introduction of Quantum Loading Scheme Ul [9|: Grover's algorithm has the func- 
tion that find the index i of a special database record record^ from the index superposition 

N-1 

of state N))- And the record record^ is the genuine answer wanted by us. However, 

i=0 

the corresponding record record^ can not be measured out unless the 1-lmapping relation- 
ship between index % and the corresponding record recordi is bound in the entangled state 

N-1 

^7=( K) \recordi)). That is, we need a unitary operation Ul to load all records that are 

i=0 

stored in a classical database into quantum state. The function of unitary operation Ul can 
be described as 

^ N-1 ^ N-1 

^=C^^ \i) |0)) \ancilla)-*—j=(y~^ j \i) \recordi)) \ancilla) (4) 



N „ VN .. 

, where ancillary state \ancilla) is known. 

Pang proposes a design method of the unitary operation Ul, that has time complexity 
0{logN) (unit time: phase transform and flipping the qubits of registers) |9|]. Operator Ul 
is so fast that its running time can be ignored when analyzing the time complexity of a 
algorithm. 

Introduction of Quantum Search Algorithm with Complex Computation (i.e., 
the Method of Rotation at Subspace) Q, H, Q, 13]: 



Grover's algorithm can find a database record according to the given index. However, 
database search is complex in general. E.g., police often hopes to find a mug shot from the 
database in which many sample photos are stored by the method of matching every sample 
photo and the photo captured by the vidicon at the entrance of airport real-time. Grover's 



algorithm is invalid for this kind of search case because the coupling between search and other 
eolation (e.g., i m a g e matching) is required at this c.e. Pang et.al. presents a quantum 



method named " rotation at subspace" [10, llll, ll2|, 113J to generalize Grover's algorithm 
to the search case with arbitrary complex computation, that is derived from the research 



of quantum image compression 13] . The method of rotation at subspace is described as 
following briefly: 

First, All input datum are stored in classical memory as database records. Assume that 
total number of records is N. All these records can be loaded into a superposition of state 
using quantum loading scheme Ul- 

Second, construct the general Grover iteration (GGI) G genera i as 

G 9eneral = (210 (ei - i)(u L y{o c yo f o c UL (5) 

, where O c denotes computation oracle such as image matching, / denotes the judge 
function (i.e., if the output of O c satisfies some conditions, let / = 1, else / = 0), Of is the 

N-1 

oracle of the judge function, and |£) = h!= K)- 

i=0 

Third, similar to Grover's algorithm, let unitary operation G genera i act on initial state 

N-1 

7)7 ( E \ i )re 9 istersl)\ )registers2 0{VN) times, We will find the Optimal Solution. 
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II. 1-DIMENSIONAL QUANTUM DFT (ID QDFT) 

First, construct the following data structure (DS) and unitary operations : 

DS1. Save DFT matrix W defined in EqfTJin classical memory as a database. And each 
line vector Wi of the matrix is regard as a record, and these records have indices 0, 1, N—l. 
DS2. Construct six registers that have data format 

| Qt)registerl \P) register2 \ % )register3 \^) register^ \ ^^i) registerb \ (Wjj • x) )register6 

. That is, 1st, 2nd, 3rd, 4th, 5th, and 6th register are used to save input parameter a, input 
parameter f3, input vector ~~x , index i, line vector Wi, and squared inner product respectively. 

DS3. Design oracle Oi nner to compute the squared inner product between vector Wi and 
x , i.e., 

|a)|0)|^|i)|w5>|O) laMl^MW^m.^) 2 ) 



We can design very simple parallel circuit to calculate squared inner product, and 
has time complexity 2t m + \log2N~\t a , where t m and t a denote the unit running time of 
multiplication and addition respectively. Because addition is more fast than multiplication 
(50t o < t m in general) and iV < 2 50 in general, the running time of Oi nner can be regarded 
as few times of multiplication (i.e., 0(l)t m ). 
DS4. Design oracle Of. 

i^> i/3> i ^f"> k> i^^;> i • -^) s > c — 1)-^<*> i^> i/5?> 1 k> i^^;> 1 c^ 7 ! • -^) s > 

, where = < 

I otherwise 

Oracle Of is used to mark the DFT coefficients Wi • ~x that satisfy the condition a < 
(Wi» ~xf < /3. 

DS5. Define ID QDFT iteration Gidqdft'- 

According to EqEl ID QDFT iteration Gidqdft is 

Gidqdft = (2|£>(£l - I){U L )-\O m ) ] O f O mner U L (6) 

Second, design the following subroutine 1 to find a coefficient Wi • ~x that 
satisfies the condition a < (W io • ~x) 2 < (3: 
subroutine 1: 

Stepl. Initialize m = 1 and set A = 6/5. (Any value of A strictly between 1 and 4/3 
would do.) 

Step2. Choose j uniformly at random among the nonnegative integers smaller than m. 
Step3. Apply j iterations of Gidqdft acting on initial state \ipo) = 

N-l 

^l«)|/?)|x)(E l<»|0>|0>. 

i=0 

Step4. Observe the 4th register: let %q be the outcome. 

Step5. Calculate value W io • ~x using classical computation. If a < (W io • ~x) 2 < (3, 
preserve iq, and exit. 

Step6. Otherwise, set m to min(\m, \^N) and go back to step 2. 

Subroutine 1 is similar to BBHT algorithm Q|, and the main difference between them is 
that Grover iteration is replaced by ID QDFT iteration Gidqdft that realizes the coupling 
between quantum search and the computation of inner product. Subroutine 1 has time 
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^ n 

condition a < (Wj • x ) 2 < (3 



complexity 0{\ j?) [8|], where M denotes the number of coefficients Wi»~x that satisfy the 



Third, design the following ID QDFT algorithm: 

N 1 

Step 0. Let AE = Wx"\\ 2 = Yl i x i x x i)i a = jf> P = ^-E, n $ = 0- We can design a 

i=0 

very simple parallel circuit to calculate value AE, and the parallel circuit has computation 
complexity 0(1) (unit time : multiplication) approximately. 

N-l 

Step 1. Generate the initial state |^o) — ~^\ a ) \P) I x )( S This can be 

i=0 

achieved in 0(log2N) steps using a |~log 2 N~\ — bit Hadamard transform, which is so fast that 
the running time can be ignored. 

Step 2. while(^ > e) 

, where e is the given threshold. 

{ _^ 

Step 2.1: Apply subroutine 1 to find a coefficient c[io] = Wi • x that satisfies the 

condition a < (c[i ]) 2 < fl] 

Step 2.2: If the DFT coefficient c[io] is the result that is not be found previously 

by this algorithm, preserve it, and let AE = AE — (c[i }) 2 , nS = nS + 1, a = j^fg , and 

(3 = AE, where a denotes the average residual energy per DFT coefficient which is not be 

obtained still by this algorithm, f3 denotes the total residual energy that is preserved by 

the DFT coefficients which are not be obtained still by this algorithm, and nS denotes the 

number of the coefficients which have been obtained by this algorithm. 

} 

It's the main idea of step 2.1 that apply subroutine 1 to find a coefficient which energy 
is bigger than the average residual energy (i.e., value a) and smaller than the total residual 

N 1 

energy (i.e., value f3). Because DFT is energy conversation transform (i.e., ^ {xi) 2 = 

i=0 

N 1 

Y (cj) 2 ), the stop criterion \r=rm > £ shows that the above algorithm will find all big 

i=0 

coefficients that preserve nearly all energy and information if threshold e is enough small. 
In addition, almost signal sequences have high redundance and the main task of DFT is 
to find and retain big coefficients to eliminate this redundance (see FigJT]). Therefore, ID 
QDFT algorithm realizes ID DFT computation approximately. 

ID QDFT has time complex O(VmN) approximately, where m denotes the number of 
the big coefficients. Because signal sequences have high redundance in general, the case 
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m <^ N often happens (see FigJT]). Thus, we regard ID QDFT has time complexity 0(y~N) 
(unit: multiplication and addition), while classical ID FFT has time complexity 0(Nlog2N) 
(unit: multiplication and addition) and even parallel DFT computation has time complexity 
O(N). QDFT can process long signal sequence at a time, while classical computer requires 
that long signal sequence must be divded into many small sections (such as 8 x 8 or 16 x 16 
section) on which DFT acts sequentially to evade the efficiency bottleneck of loading data 



III. 2D QDFT 

2-D DFT (see Eqj3]) is a separable linear transformation. That is, the result of 2-D DFT 
may be obtained by first taking transforms along the columns of F and then along the rows 
of that result |l|. That is, C = W N FW N = (W N F)W N . We define 

G = W N F = (Wi • Tj)nxn 

, where Wi is the line vector of matrix Wn and fj is the column vector of input matrix F. 
The main task of 2D QDFT is to calculate out matrix G and matrix GWn (i-e., C). 2D 

QDT is described as following: 

First, Construct the following data structures (DS) and unitary operations: 

DS1: Store all of line vector Wi in a database, and each index of line vector is denoted 

by i. Store all of column vector fj of matrix F in the database, and each index of column 

vector is denoted by j. 

DS2: Construct seven registers that has data format |a)|/3)|2)|j)|H / j)|/j)|(W / i • fj) 2 }- 
DS3: Design two loading operation Uli and Ul2 according to Eqfj] and ref. where 

Uli and Ul2 are applied to load vector Wi and fj into registers from classical database 

respectively. That is, 

\a}\m\m\m u ^\a}\m\mwm\o) 

l«)|/5)K)|j)|0)|0)|0)^ 2 |«)|/5)|z)|j)|0)|^)|0) 
DS4: Design oracle Bi nner to calculate squared inner product. That is, 



l«> \P) l*> \JWM} |0> \a) \(3) |i> \j)\W)\7 3 ) I (Wi . £) 2 > 



DS5: Design oracle O'j : 

\a)\m\3)\tfi)\7j)0i • tf) 5 (-i)^i«>i^iou>iw5>i^>i(w5«^) 2 > 

, where = 



l*/a<(WW,) 2 </3 

otherwise 
DS6: Design 2D QDFT iteration Gidqdft'- 



g 2DQDFT = (2irxri - ^(t/Li)^^)^^^)^/^^^^! 

iV-1 AT-1 

, where |£'> = £ £ £ IW^- 

N-1 N-1 

Second, Generate the initial state |^ ) = jj £ £ \ a )\P) l )l )l )- And a PP!y 
the method of ID QDFT, we will calculate out all elements of matrix G. Apply the same 
method again, we can calculate out all elements of matrix C. 

2D QDFT has time complexity O(N), while classical 2D FFT has time complexity 
0(N 2 log 2 N). 



IV. USING ID QDFT TO CALCULATE CONVOLUTION 

For the given two N-dimensional vectors it = (u , u±, tt/v-i) and v = (v , v±, Uzv-i), 

N-1 

let Wk = £ u j v (k-j)modN, where < k < N. Then vector w = (w ,wi, ...,Wjv-i) is called 
the periodic convolution of vectors u and v . We often denote the convolution by 
symbol ~w = u * v . If ~u and have different dimensions, add zero components to each 
vector such that both vectors have the same dimension. The periodic convolution of the 
new vectors can be regarded as the result. Algebraically, convolution is the same operation 
as multiplying the polynomials whose coefficients are the elements of ~u and ~v 

Holding the view of signal processing, we can regard vector u as input signal sequence 
and vector is the function of a physical apparatus. After vector u passing through 
the apparatus, convolution w is came out. Furthermore, the output sequence Tu becomes 
smoother than the input sequence ~u in general. Therefore, convolution is very important 
for signal or information processing 



The convolution theorem [l] says, roughly, that convolving two sequences is the same as 
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multiplying their Fourier transforms. That is |l| 



DFT(~u * ~v) = DFT{~u)DFT{~v) (7) 

In general, the big coefficients of DFT(v) are few and other coefficients have little 
contribution to the result. That is, we can calculate only big coefficients of DFT(u) 
and DFT(v) and discard small coefficients to speedup the computation of convolution. 
According to this property, quantum algorithm for convolution estimation is described 
as following: 

Stepl. Calculate DFT(Tl) and DFT(~v) using the ID QDFT. 
Step2. Calculate DFT(~u * v) according to Eqd 

Step3. Calculate the inverse of DFT(~u * v) to obtain convolution ~u * v . 

It should be noticed that the efficiency of the above algorithm depends on the redundancy 
(or smooth property) of the signal sequences. The more higher the redundancy is, the more 
fast the quantum algorithm runs. Fortunately, many digital signal sequences have high 
redundancy, that is just the existence reason of modern signal processing technique. 
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